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Abstract Gravitational waves are perturbations in the 
spacetime that propagate at the speed of light. The 
study of such phenomenon is interesting because many 
cosmological processes and astrophysical objects, such 
as binary systems, are potential sources of gravitational 
radiation and can have their emissions detected in the 
near future by the next generation of interferometric 
detectors. Concerning the astrophysical objects, an in¬ 
teresting case is when there are several sources emitting 
in such a way that there is a superposition of signals, re¬ 
sulting in a smooth spectrum which spans a wide range 
of frequencies, the so-called stochastic background. In 
this paper, we are concerned with the stochastic back¬ 
grounds generated by compact binaries (i.e. binary sys¬ 
tems formed by neutron stars and black holes) in the 
coalescing phase. In particular, we obtain such back¬ 
grounds by employing a new method developed in our 
previous studies. 

Keywords Gravitational waves • Stochastic back¬ 
ground • Coalescing binaries • Compact objects 

1 Introduction 

Binary systems are among the best known sources of 
gravitational waves. Specially interesting are the dou¬ 
ble neutron star systems (NSNS systems), the binaries 
formed by a neutron star and a black hole (BHNS sys¬ 
tems) and the binaries of black holes (BHBH systems), 
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because their emissions have a high probability of being 
detected in the near future. 

Besides, it is well known that binaries lose energy 
and momentum via emission of gravitational waves, 
which cause reductions in their orbital distances, con¬ 
sequently increasing the orbital frequencies. Moreover, 
for systems in circular orbits, the frequency of the emit¬ 
ted waves is twice the orbital frequency. The process 
continues until the systems reach the coalescing phase, 
where the systems leave the periodic regime and start 
the merging phase. 

On the other hand, concerning the study of the grav¬ 
itational radiation itself, the stochastic backgrounds are 
of special interest. Backgrounds can be generated when, 
for example, there is a superposition of signals of several 
sources, resulting in smooth-shaped spectra spanning a 
wide range of frequencies. In particular, we are con¬ 
cerned in this paper with the backgrounds generated 
by a population of coalescing compact binaries formed 
from redshifts ranging from z = 0 up to z ^ 20, i.e., 
cosmological binaries. 

We can find in the literature some very interest¬ 
ing works on this issue, such as [TJ[5J[31I1] , where the 
authors calculated the backgrounds generated by coa¬ 
lescing NSNS systems. Generally speaking, they used 
Monte Carlo techniques to simulate the extragalactic 
population of compact binaries. Besides, these authors 
considered the time evolution of the orbital frequency 
by using the “delay time” (that is, the interval of time 
between the formation and the coalescence of the sys¬ 
tems) in their calculations. In our method, we consider 
the time evolution in an explicit form, where the equa¬ 
tion describing the evolution of the orbital frequency is 
taken at each instant of life of a given binary. 

In Zhu et al.[S], the authors consider the spectra 
generated by coalescing BHBH systems. In this paper. 
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they assume average quantities for the energy emissions 
of single sources; on the other hand, in our calculations 
we consider all the values for the orbital parameters a 
system can have, in the form of distribution functions. 

Still concerning the spectra generated by coalesc¬ 
ing BHBH systems, Marassi et al. adopted an up¬ 
dated version of the SeBr0 population synthesis code, in 
which the masses of the black holes range from ~ GMq 
and ~ 20Mq. 

Roughly speaking, the various papers cited above 
show results that are characterized by spectra with fre¬ 
quencies ranging from ~ lOHz and ~ lO^Hz and with 
maximum amplitudes located in the interval ranging 
from ^ 400Hz to ^ 800Hz. As shown below, the back¬ 
grounds we generated have similar forms to the ones we 
mentioned, though our results show, in general, higher 
amplitudes. This difference will be discussed timely. 

Further, as it will be shown, our method has the very 
useful characteristic of being numerically simple, in the 
sense that it does not demand a heavy computational 
work. 

In this paper, the spectra generated by the coalesc¬ 
ing binaries will be calculated by means oflZllH] 

^BG = “ [ ^sourcedR, (1) 

^obs J 

where /ibg represents the dimensionless amplitude of 
the spectrum, z/obs is the observed frequency, hgource is 
the amplitude of the signals generated by each source 
and dR is the differential rate of generation of gravi¬ 
tational waves. It is worth pointing out that o was 
deduced from an energy-flux relation. In fact, in a pa¬ 
per by de Araujo et al.[5], the authors gave a detailed 
derivation of this equation, showing its robustness. Ac¬ 
tually, one can use m in the calculation of differ¬ 
ent types of stochastic backgrounds, provided that one 
knows ^source and the corresponding dR to the case one 
is dealing with. 

Here, ^source has the formPITU] 

^source = 7.6 X 10“^^ 



where v is the emitted frequency that, in this case, is 
the frequency emitted by a coalescing system (note that 
ly and z/obs are related to each other by ini), fi is the 
reduced mass and M is the total mass. The differential 
rate dR in writing in following form 

= ^^h{'^i)f2{rn2)dzdmidrn2, (3) 

dV dz 

^ see WWW.sns.ias.edu/~ starlab 


where dV is the comoving volume element, ^ is the red- 
shift and /i(toi) and / 2 (to 2 ) are the mass distribution 
functions of the components of the systems. In fact, for 
neutron stars such distributions are given by Dirac’s 
delta functions, since we are considering that all these 
objects have the same mass of 1.4 Mq; for black holes 
we use the function [TT] 

/(m) = 0.332exp [-0.347(m - 7.8)^] , (4) 

where m is given in solar mass units. 

The term dV/dz is known from Cosmology, and the 
problem of determining dR comes down to the calcula¬ 
tion of dR/dV. 

At this point, one could ask in what the present 
study differs from the previous ones. We will see that 
the main difference has to do with the application of a 
new method to calculate dR developed in our previous 

papers [T211I3]- 

The paper is organized as follows: in Section ^ we 
show the main steps to obtain dR; then, with ([T|), dR 
and at hand, we calculate /ibg for the three families 
of compact binaries, which are explained in Section[3] in 
Section m we present the results and discuss, in partic¬ 
ular, the detectability of the backgrounds studied here 
by the interferometric detectors Laser Interferometer 
Space Antenna (LISA) (now evolved LISA (eLISA)), 
Big Bang Observer (BBO), DECI-Hertz Interferome¬ 
ter Gravitational wave Observatory (DECIGO), Ad¬ 
vanced Laser Interferometer Gravitational wave Ob¬ 
servatory (ALIGO), Einstein Telescope (ET), and the 
cross-correlation of pairs of ALIGOs and ET; finally, in 
Section Owe present our conclusions. 


2 Calculation of dRfdV 

We present here the main steps for the calculation of 
dR/dV (we refer the reader to Refs. [T^fT^ for details). 
We start by writing dR/dV in the form 

dR dv f 

dy^H «binH(i/,t,to)ato, (5) 

where the expression in the right-hand side is the for¬ 
mation rate of systems per comoving volume that reach 
the frequency z/ at the instant t (see the derivation in 
Appendix E] and also in Ref. m)- Also, to refers to the 
instant of birth of the systems. 

In (O, iL(z/, t, to) is the frequency distribution func¬ 
tion of the binaries, which has the form (we refer the 
reader to Appendix [B] and Ref. 1 12) for the derivation of 
this equation): 
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where z/q is the initial frequency; C, cr and f are con¬ 
stants given in Table [2] and r is the orbital distance, 
related to v by means of Kepler’s third law. 

On the other hand, ribin is the binary formation rate 
density that, for NSNS systems is given bys 

P*{Zo) -liiT -3 

rinsns = Ansns-r-- yi MpC . (7) 

1 + 2^0 

Here, zq refers to the redshift of birth of the systems, 
which is related to to via the usual expression that can 
be found in any textbook of cosmology; and p*{zo) is 
the star formation rate density (SFRD). 

There are, in the literature, many different proposals 
to the SFRD, although they do not differ from each 
other very significantly. Here we adopt, as a fiducial 
one, that given by Springel and Hernquist [2] , namely 

P*{z) = pm-^—^ . , MQyr-^Mpc~^ (8) 

where a = 3/5, /3 = 14/15, Zm = 5.4 and with pm = 
O.15M0yr“^Mpc“^ fixing the normalization. 

Besides, in o ; -^nsns — Pn.s fp<Pns (see m) is the mass 
fraction of stars that is converted into neutron stars, 
where ,5ns is the fraction of binaries that survive to the 
second supernova event; fp gives the fraction of massive 
binaries (that is, those systems where both components 
can generate supernovae) and ^ns is the mass fraction 
of progenitors that originates neutron stars, which, in 
the present case, is calculated by 

^25 

= / 4>{m)dm (9) 

Js 

where />(m) = is the Salpeter mass 

distribution [IS] with x = 1.35 and A = 0.17. Numer¬ 
ically, we have /3ns = 0.024, fp = 0.136 and ^ns = 
5.97 X IO-^Mq^ 

Now, considering © , we can write ([S]) in the form 


dR _ 


nsns 


P*,c^z), 


( 10 ) 


where, following the notation adopted by Zhu et alia, 
p*^c{z) is calculated by 


P*Az) 


P*{zo) 

1 + 2:0 


. , dv 


dta. 


( 11 ) 


In this context, following Wu et aim, Ansns is given 
by Ansns = ?'o,nsns/p*,c(0) where ro_nsns is the local coa¬ 
lescence rate and /9*,c(0) has the form 


P*,c(0) 


P*{zo) 

1 + 2:0 




dv 

dt 


dtn 


z=0 


( 12 ) 


For BHNS and BHBH systems, we have similar ex¬ 
pressions, but with different values of noin- Such values 
can be estimated by means of the results found in Bel- 
czynski et al. m. where these authors claim that the 
population of binaries is formed by 61% of NSNS, 30% 
of BHBH and 9% of BHNS binaries. Therefore, nohns 
and nbhbh can be related to rinsns by means of these 
proportions. 

It is worth mentioning that Belczynski et al. m 
studied compact binaries with merger times lower than 
IQiOyr; jg^ consldered coalescing binaries. Ba¬ 
sically, in the simulations they used, the binaries are 
formed through several different channels. Specifically, 
NSNS systems are formed through 14 different chan¬ 
nels, where there is a predominance of channels contain¬ 
ing hypercritical accretion between a low-mass helium 
giant and its companion neutron star. On the other 
hand, BHNS and BHBH systems are formed through 
just for and three channels, respectively, where there is 
a moderate predominance of mass transfer events. 


3 Calculation of the spectra 

In this Section, we present the calculations of the spec¬ 
tra generated by NSNS, BHBH and BHNS systems. Al¬ 
though we are using m for the three cases, the calcula¬ 
tions are different for each family of binaries. Therefore, 
we will show the calculations separately. 


3.1 NSNS systems 

Neutron stars, according to theories of stellar evolution 
and observations (see, e.g., Ref. m) have characteristic 
masses that fall in a narrow interval around 1.4M0. So, 
in this paper we are considering that all neutron stars 
have masses of 1.4M0, which is a realistic choice, and 
at the same time, a simplification. 

One could ask how the results would be affected by 
the choice of the neutron star (NS) equation of state 
(EOS). Since we are considering the coalescing phase, 
the results depend mainly on the mass of the NSs. 
NSNS systems are characterized by a specihc coales¬ 
cence frequency that, according to, e.g.. Ref. [IH], may 
be considered as « 900Hz for a pair of two 1.4M0 NSs. 
The choice of the EOS is certainly important for the 
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subsequent phase of evolution of the system, when the 
merger phase takes place. 

Although there is just a value for the coalescence fre¬ 
quency, the spectrum will be spread over a wide range 
of frequencies. This behavior is due to the cosmic red- 
shift, because systems emitting at the same frequency, 
but at different redshifts, will generate signals with dif¬ 
ferent observed frequencies, obeying 

^obs — z ; (1^) 

1 + z 

where, in this case, we are considering ly = 900Hz. 
Moreover, as we are considering that the redshift has 
the minimum and maximum values given by ^min = 0 
and Zrnax = 20, respectively, the observed frequencies 
will have minimum and maximum values of = 

42.86Hz and = 900Hz, respectively. 

Since the masses of the components of the systems 
and their coalescing frequencies are the same, the spec¬ 
tra will only depend on the redshift 2 . However, from 
m one notices that there will be just one value of r'obs 
for each value of 2 , such that © must be handled in a 
particular way. First, let us rewrite o in the form 

2 ^ j-Zp+Sz 

^BG = “ / f^source^R (14) 

Vohs J Zp — Sz 

where 5z << 1, and Zp is the value of redshift corre¬ 
sponding to each observed frequency by means of (HSl). 

Now, in order to calculate o, dR can be written 
in the following form 

dR= ^V{zp)S{z - Zp)dz (15) 

where V(zp) is the value of the comoving volume at Zp 
and S{z — Zp) is the Dirac’s delta function. Therefore, 
substituting (fTSl) in m and integrating, we have 


where risco is the radius of the ISCO of the black hole, 
which is related to the Schwarzschild radius, namely 


Dsco — drgchw 


6GTObh 


Substituting (ITS)) in (IT71) we have 


(18) 


tGV^' 


rribh 


m- 


bh 


(19) 


Since we are considering that black holes have masses 
in the interval SMq — 2 OM 0 , the emitted frequen¬ 
cies will have minimum and maximum values given by 
i^min _ ]^x4Hz and = 997Hz, respectively. More¬ 
over, considering (dsi) and (di), one notes that for each 
value of Vohs, one has continuous ranges of values for z 
and TObh- 

However, in the case of BHNS systems, m must be 
integrated over z but not over rrihh, since these variables 
are not independent. Moreover, for a given value of t'obs 
one needs to determine the maximum and minimum 
limits of integration, which are given by 

^min 

2min =-1 (20) 

^obs 

^max 

— 1- (^1) 

^obs 

When ^min < 0 or 2:max > 20, we set Zmin = 0 and 
Zmax = 20, respectively. Once the interval of integration 
is set, the variables mbh a-nd u in the integral can be 
written as functions of z by solving m and m- 
In addition, note that in this case one needs to con¬ 
sider in the calculation of the distribution function 
given by As a result, we have 

I y-zmax dRdV , , , , 

^BG = — / Kource^^flMdZ- (22) 

^^obs ./Zjjjjjj dz 


h'^ - 

^BG — 


l^ohs 


, n dR 

-h^ - 


Vizp). 


Zp,l^obs 


(16) 


3.2 BHNS systems 


The study of the coalescence of BHNS systems is more 
complicated than the case of NSNS systems, because 
the frequency of coalescence depends on the mass of 
the black hole. It is usually assumed that the coales¬ 
cence occurs when the neutron star reaches the inner¬ 
most stable circular orbit (ISCO) of the black hole. So, 
using Kepler’s third law and recalling that the emitted 
frequency is twice the orbital frequency, we have 



G(mns + TObh) 


(17) 


3.3 BHBH systems 

In this case, the frequency of coalescence depends on the 
values of the two components of the system. Following 
Marassi et al.[5], i/ is given by 

+ boT] + Co 

■'=0 -iSJ- 

where 77 = TObhi TObhj/-AT^ is the symmetric mass ratio 
and M is the total mass. The polynomial coefhcients are 
flo = 2.974 X 10“^, bo = 4.481 x 10“^ and cq = 9.556 x 
10“^. According to (1^ . for each value of v, there will 
be a continuous set of pairs of values for TObhi and TObh 2 
that satisfies this equation. Moreover, the masses are 
not independent, in fact they are related to each other 
by (1^51) . Therefore, we should integrate over one of the 
masses (or one parameter which describes both masses.) 
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However, it is not possible to solve (ESI) analytically, 
i.e, to write the masses as functions of v. Therefore, we 
need to use approximations. First, consider that mbhi 
and mbh 2 (which we refer as to mi and m 2 ) are related 
to each other by 

mi = km2 (24) 


where the variable k is greater than or equal to one. 
Substituting (IM)) in (OSl) we can write the masses mi 
and m 2 as functions of v and k, namely 


mi = k 


£_J_ 

G TTU 
:3 1 


Oo 




bnk 


m 2 = 


G 


TTiy 


(1 + k)^ 

aok^ 


(1 + k)^ 

bok 


CQ 

1 + k 
Co 


(1 + fc )5 (1 + fc)® \ + k 


A suitable approximation for (I2dl) is given by 


mi + m 2 


c^ Co 


G 


(25) 


since 77 < 1 for all values of mi and m 2 . Considering 
that m 2 = minin = 5 hl 0 , we use (1251) in order to es¬ 
timate a hrst value for mi; next, we take this pair of 
values for the masses and calculate a first value for ry; 
we then correct the value for mi using 

0077^ + b^T] + Co , , 

mi-|-m 2 Ri—-. (26) 

Lr TTV 

With this new value of mi, we repeat the above process: 
we calculate rj again and use (1261) , bearing in mind that 
this process may be performed an arbitrary number of 
times in order to yield more accurate values for mi. In 
the cases where we have mi > mmax = 2 OM 0 at the 
end of the process, we consider mi = mmax and perform 
an analogous process to find out m 2 . 

With the pair (mi,m 2 ) at hand, we calculate the 
maximum value of k: 
h _ ^ 

^max — 

7712 

Thus, in order to cover all possible values of mi and 
m 2 , we consider that fc, in (I25L ranges form one to 
k 

'^max- 

Finally, (HD can be written as follow 


(27) 


^bg — 


1 /*fcmax /*2max 

— / / 

I I "source 

^obs h 

X^^fii'mi)f2im2)dzdk. 
dV dz 


(28) 


4 Results and discussion 

As usual, in the literature, we represent the back¬ 
grounds in terms of the strain amplitude Sh, which is 
given by 

Sn = ^ (29) 



Fig. 1 Backgrounds for NSNS, BHNS and BHBH systems 
in terms of the strain amplitude 



Fig. 2 Backgrounds for NSNS, BHNS and BHBH systems 
in terms of the energy density parameter 


and also in terms of the energy density parameter I7gw, 
which reads (see, e.g.. Ref. [19]) 


12 ™ = 


ml 


BG- 


(30) 


The spectra are shown in Fig. [T] and Fig. [5] 
for Sh and 12gw, respectively. Note that these spec¬ 
tra are mainly compared to the sensitivity curves 
of ET|20] and ALIGO|20], since the spectra are 
in their frequency bands. The sensitivity curves of 
the proposed space-based antennas LlS^jj, eLISApi]. 
BBO[52], DECIGO|53] are also shown, but given their 
frequency bands they cannot detect the spectra studied 
here. 

From Fig.jT] one notices that the background gener¬ 
ated by the three families of compact binaries are below 
the sensitivity curves of the interferometric detectors. 
Besides, one notices that the background generated by 
the BHNS systems have higher amplitudes when com¬ 
pared to the ones generated by NSNS and BHBH sys¬ 
tems. On the other hand, the spectrum corresponding 
to BHBH systems has the lowest amplitudes. 


^obs 


see WWW. srl.caltech.edu/"shane/sensitivity 
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Since our calculations depend on some parameters 
and functions, it is worth investigating how our results 
are affected by different choices of these quantities. 

First, let us consider the masses of the components: 
in ([2), if we multiply both masses by a factor of q, ^source 
will be multiplied by a factor of ~ Important 

variations in the amplitudes would occur only if g << 1 
or g >> 1. 

From o, m and (1^51) . one notices that \/Sh oc 
\/ P*{zo) and y/Sh oc VAnsns- Therefore, if ones multi¬ 
ply p*{zo) or Ansns by, say, a factor of 10, the amplitudes 
shown in Fig.[T]will increase by a factor of ~ 3.2. There¬ 
fore, for realistic scenarios, different choices for p»(2o) 
and Ansns would have small effects on the amplitudes of 
the backgrounds. 

Comparing Fig. [5] with similar studies found in the 
literature (see, e.g., Refs. [Il[2j|3l|4j[5l[6]), one sees a good 
agreement concerning the shapes of the spectra, al¬ 
though our results show higher amplitudes. For ex¬ 
ample, in m one sees that ^ 10 at 

i^ohs ~ 800Hz for the backgrounds generated by NSNS 
systems, while for our corresponding spectrum, shown 
in Fig. [U we have y/Sh ^ 5 x 10“^®Hz~^/^ at the same 
frequency. 

Comparing our results for BHBH systems (see 
Fig. [2]) with the results found in , one can note some 
similarities: the amplitudes increase until a maximum 
value in the range 400 — 500Hz and then they have a 
sharp decrease. Concerning the amplitudes, we have a 
maximum value of ^ 5 x 10“®, while in Zhu et al the 
value is ~ 2 X 10“®. 

Marassi et al [5] also study backgrounds generated 
by BHBH binaries. These authors discussed different 
models, and the resulting spectra present maximum 
amplitudes ranging from 10’ < f2gw < 5 X 10 ® for a 

frequency band around ^ 500Hz. 

It is worth mentioning that, generally speaking, the 
spectra is model dependent. Therefore, different as¬ 
sumptions lead to different backgrounds. In Ref.[T], for 
example, the population of binaries is such that the 
maximum probability of coalescence is around z = 1.4. 
Therefore, for ^ < 1.4 there is a relatively small propor¬ 
tion of coalescing systems emitting; in our calculations 
we do not consider such a behavior. This difference in 
the proportion of systems at lower redshifts could ex¬ 
plain our higher amplitudes as compared to the ones of 
Ref. [T]. 


4.1 Cross-correlation of pairs of detectors 

Although the spectra (signals) shown in Fig. [1] are be¬ 
low the sensitivity curves of the detectors, it could well 


Table 1 S/N for the three families of coalescing compact 
binaries, considering T = lyear 


System 

ALIGO 

ET 

NSNS 

6.3 X 10-3 

5.0 X 102 

BHNS 

3.4 X 10-2 

1.2 X 103 

BHBH 

1.6 X 10-3 

6.4 X lO^ 


be possible detect them by correlating the outputs of 
two or more detectores. For the correlation of two in¬ 
terferometers, the detectability of a given signal can be 
quantified by means of the so called signal-to-noise ratio 
(S/N), namely PSII^T] : 


[S/Nf 





(31) 


where S\ and 5"^ are the spectral noise densities, T is 
the integration time, and ^{y) is the overlap reduction 
function, which depends on the relative positions, spa¬ 
tial orientation, and distances of the detectors; and I7gw 
is given by o. 

In Table [T] one can see the S/N for the three families 
of compact binaries, in particular for pairs of ALIGOs 
and ET. 

From Table[l] one notices that ET could in principle 
detect the backgrounds where the spectrum generated 
by BHNS systems would have higher probability of de¬ 
tection; for pairs of ALIGOs, the low values of the S/N 
ratio indicate a non detection. 


5 Conclusions 

In this paper, we calculate the stochastic background 
of gravitational waves generated by coalescing compact 
binaries, using a new method developed in our previous 

studies mm- 

We show that, of the three spectra considered in 
this paper, the one generated by BHNS systems has 
the highest amplitudes, while the background by BHBH 
systems show the lowest amplitudes. Moreover, one no¬ 
tices slight differences in the forms of the spectra, which 
are due to the different methods used to calculate them 
by means of O- 

We found that the backgrounds calculated here 
would not be detected by the interferometric detectors 
such as LIGO and ET, although thanks to the cross¬ 
correlation of signals ET could, in principle, detect such 
signals. Particularly, we found that the spectrum gen¬ 
erated by BHNS systems have the highest S/N ratio, 
while the one corresponding to BHBH systems presents 
the lowest S/N. 

Concerning the dependence of our results on the 
parameters used in the calculations, we found that the 
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masses of the components of the binaries, as well as 
P*{zq) and Ansns, do not strongly influence the back¬ 
grounds. Besides, a particular choice for the NS EOS 
does not affect the results either. 

We compared the spectra studied here with some 
interesting results found in the literature. One notices 
similarities in their shapes, namely: maximum frequen¬ 
cies of ^ lO^Hz and maximum amplitudes in the range 
400 — 800Hz. Roughly speaking, these characteristics 
are common to the three families of binaries. 

On the other hand, our amplitudes given in terms 
of \fSh are in general higher than the ones found in 
the literature by roughly one order of magnitude. We 
concluded that such a difference is mainly due to pop¬ 
ulation characteristics assumed. Therefore, generally 
speaking, the spectra is model dependents. 
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A Appendix 


In this appendix, we present the main steps of the derivation 
of dRjdV. For further details, we refer the reader to Refs. 

Tmm. 

We derive dRjdV by means of an analogy with a prob¬ 
lem of Statistical Mechanics. In this problem, the aim is to 
calculate the number of particles that reach a given area A 
in a time interval dt, i.e., the objective is to calculate the flux 
F of particles. Basically, this flux is calculated by counting 
the particles inside the volume dV = Adx, adjacent to the 
area A, that are moving towards A with velocity v = dxjdt, 
where v obeys a distribution function p{v). Hence, the flux is 
obtained by integrating over all the positive values of v. 

With some modifications, this method can be used to de¬ 
termine dR/dV. First, we substituted the spatial coordinate x 
by the frequency u and the velocity v by the time variation of 
the frequency, which is defined by = dvJdt. Therefore, the 
number of systems in the interval dv adjacent to a particular 
frequency v is given by 




f ifi{v')du' 


(32) 


where ip{v) is the non-normalized distribution of frequencies. 

Considering that the distribution rj[vv) gives the number 
of systems which have Vi, in the interval dv^, the number of 
systems in dv and with values of Vi, in the interval dvu is 
given by 


Now, the next step is to determine ip{v) and ri{vv). First, the 
distribution ^p(v) is written in the form 

ipiv) = J n^ij^{to)H(v)dto (34) 


where to is the instant of birth of the systems, is the 

formation rate density of the NSNS, and H{v) is given by 

( 0 . 

In the derivation of (I34II . we consider initially H{v), from 
which we have 


dn = H(v)dv, 


(35) 


which is the fraction of systems originated at the time to and 
that have frequencies in the interval dv. Now, using we 

can write explicitly 


dn 

Now, integrating over dto, we get 


dn 

dV 


j nbin(fo)H(iz)dfo 


dv. 


(36) 


(37) 


where the expression in brackets is the number of systems 
per unit frequency interval and per comoving volume, which 
is the desired distribution function '.p{v). 

On the other hand, rj{v^) will have a peculiar form. First, 
note that the derivation of iliil) yields 

dv 11 

Vv = — ccv 3 (38) 

dt 

after some algebraic manipulations. We conclude that there 
will be just one value of Vi, for each value of v, which allows 
us to write rj{v^) as a Dirac’s delta function, namely 


ri{v^) = NS{vv — v^^p), (39) 

where N is the total number of systems and is the par¬ 
ticular value of Vi, corresponding to each frequency v. 

Notice that the denominator of the term between paren¬ 
thesis in 13211 is the total number of systems. Now, using the 
function given by iP9)) and changing the differential dv by 
means of the chain rule, 13311 assumes the form 

/ ip(v)^dt\ 

dp = j -—- NS{v^ — v^^p)dv^. (40) 


Integrating over v^, and rearranging the terms, we obtain 


„ , .dv 

R = ^{u)- 


(41) 


where R is the number of systems per time interval dt. Re¬ 
calling that the rate R is per comoving volume, one has 


dv dR 
' dt dV 


(42) 


B Appendix 

In this paper, the frequency distribution was derived from the 
semi-major axis distribution given by the following gaussian 
function [TB] 


/g(^) = C'exp 


— (r — r)^ 


2o-2 


(43) 


where r is the semi-major axis and the parameters r, C and 
cj have the values given in Table [2] 
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Table 2 Parameters of the distribution functions of the or¬ 
bital distances 


System 

C 

f(R0) 

<t{Rq) 

NSNS 

0.070 

0.6 

0.2 

BHNS 

0.015 

5.5 

1.5 

BHBH 

0.070 

11 

2.5 


First, we changed variables via f{r)dr = g{Q)dQ with 
the aid of Kepler’s third law, where Q is the angular orbital 
frequency. Note that Q depends on time (see Ref. m), namely 


n = 





- 3/8 


(44) 


where K = 96mim2 5c®G®/®(mi -|- rui and m 2 

are the masses of the components of the system and flo is 
the initial frequency. So, carrying out a change of variables 
via g{Qo)dQo = H{Q)dQ, where 17o was associated with the 
variable Q in g{0), one has 


H{Q) = 




'GM' 

Att 


1/3 


— (r — r)^ 

2 ^^ 


(45) 


17. B. W. Carroll and D. A. Ostlie, An Introduction to Mod¬ 
ern Astrophysics (Addison-Wesley,2007) pp 569, 578, 639 

18. G. Poghosyan, R. Oechslin, K. Uryu and F K Thiele¬ 
mann, Mon. Not. R. Astron. Soc. 349, 1469-80 (2004) 

19. V. Ferrari, S. Matarrese and R. Schneider, Mon. Not. R. 
Astron. Soc. 303, 258 (1999) 

20. C. K. Mishra, K. G. Arun, B. R. Iyer and B. S. 
Sathyaprakash, arXiv:1005.0304v2 (2010) 

21. P. Amaro-Soane et al, arXiv:1202.0839v2 (2012) 

22. C. Cutler and J. Harms, Phys. Rev. D 73, 042001 (2006) 

23. K. Yagi and T. Tanaka, arXiv:0908.3283v2 (2010) 

24. J. M. Lattimer, Annu. Rev. Nucl. Part. Sci. 62, 485 

( 2012 ) 

25. B. Kiziltan, A. Kottas, M. De Yoreo and S. E. Thorsett, 
Astrophys. J. 778, 66 (2013) 

26. B. Allen, Relativistic Gravitation and Gravitational Ra¬ 
diation (Cambridge Univ. Press, Cambridge,1997) 

27. B. Allen and J. D. Romano, Phys. Rev. D 59, 102001 
(1999) 

28. P. C. Peters, Phys. Rev. 136, B1224 (1964) 


Finally, it would be necessary to perform a further coordinate 
transformation in order to write H{f2) as a function of the 
emitted frequency v. Such a transformation, calculated by 
means of 1? = ttu, is trivial. Besides, r is written as a function 
of Q by means of Kepler’s third law. 
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